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According to general relativity a perturbed black hole will settle to a stationary configuration by the emis- 
sion of gravitational radiation. Such a perturbation will occur, for example, in the coalescence of a black hole 
binary, following their inspiral and subsequent merger. At late times the waveform is a superposition of quasi- 
normal modes, which we refer to as the ringdown. The dominant mode is expected to be the fundamental mode, 
I = m = 2. Since this is a well-known waveform, matched filtering can be implemented to search for this signal 
using LIGO data. We present a search for gravitational waves from black hole ringdowns in the fourth LIGO 
science run S4, during which LIGO was sensitive to the dominant mode of perturbed black holes with masses 
in the range of 10A/ Q to 500 M @ , the regime of intermediate-mass black holes, to distances up to 300 Mpc. 
We present a search for gravitational waves from black hole ringdowns using data from S4. No gravitational 
wave candidates were found; we place a 90%-confidence upper limit on the rate of ringdowns from black holes 
with mass between 85 Mq and 390 Mq in the local universe, assuming a uniform distribution of sources, of 
3.2 x 10~ 5 yr _1 Mpc~ 3 = 1.6 x 10" 3 yr" 1 L^, where L 10 is 10 10 times the solar blue-light luminosity. 



PACS numbers: 04.30.-w, 04.30.Tv, 97.60.Lf, 04.80.Nn 

I. INTRODUCTION 

The existence of intermediate mass black holes, IMBHs, 
(2OM < M < 10 6 M o ) has been under debate for several 
decades. While general relativity does not preclude IMBHs, 
there had been no observational evidence for their existence 
until recently. Electromagnetic observations have indicated 
that ultra-luminous X-ray sources, that is, sources radiating 
above the Eddington luminosity for a stellar mass black hole, 
may be powered by IMBHs. Strong evidence in support of 
this argument has recently been reported with the discovery 
of a source whose luminosity implies the presence of black 
hole with mass of at least 500 Mq Q~). Further information 
pertaining to IMBHs may be found in recent comprehensive 
review articles fl2j[3]. 

Predictions have been made for the rate of detection of 
ringdowns from IMBHs in Advanced LIGO. Ref. H pre- 
dicts a rate of 10 events per year from IMBH-IMBH binary 
coalescences. When scaled to the sensitivity of the data set 
under consideration in this investigation, this prediction be- 
comes 10~ 4 yr _1 [5]. Ringdowns following coalescences of 
stellar-mass BHs with IMBHs could also be detectable with 
Advanced LIGO, with possible rates of tens of events per year 

Detection of gravitational radiation from IMBHs however, 
would provide unambiguous evidence of their existence. In 
order for such an object to reveal itself through gravitational 
radiation it must come to be in a perturbed state, for exam- 
ple as the remnant of the coalescence of two IMBHs. Current 
ground-based gravitational wave detectors, such as the Laser 
Interferometer Gravitational- Wave Observatory (LIGO), op- 
erate in an optimal frequency range for the detection of the 
ringdown phase of the binary coalescence of IMBH binaries. 
In this paper we describe a search for ringdown waveforms in 
data from the fourth LIGO science run, S4. 



II. THE RINGDOWN WAVEFORM 

A series of studies within a linearized approximation to Ein- 
stein's equations and also full-blown numerical simulations 
have shown that the gravitational wave signal from a perturbed 
black hole consists of roughly three stages Q: (i) A prompt 
response at early times, which depends strongly on the initial 
conditions, and is the counterpart to light-cone propagation; 
(ii) An exponentially decaying "ringdown" phase at interme- 
diate times, where quasi-normal modes, QNMs, dominate the 
signal, which depends entirely on the final black hole's pa- 
rameters; (iii) A late-time tail, usually a power-law fall-off 
of the field. The ringdown phase, which is the focus of this 
work, starts roughly when the perturbing source reaches the 
peak of the potential barrier around the black hole, and con- 
sists of a super-position of quasinormal modes. For instance, 
during the merger of two black holes, the start of the ring- 
down is roughly associated with the formation of a common 
apparent horizon, which also corresponds to the peak of the 
gravitational-wave amplitude. For black holes in the LIGO 
band this is on the order of tens of milliseconds after the in- 
nermost stable circular orbit of the binary. Each quasi-normal 
mode has a characteristic complex angular frequency u>i m ; the 
real part is the angular frequency and the imaginary part is the 
inverse of the damping time t. Numerical simulations (for ex- 
ample [8 1) have demonstrated that the dominant mode is the 
fundamental mode, I = m = 2, and that far from the source 
the waveform can be approximated by 

ho^^n^A^e-^ (1) 

where A is the dimensionless amplitude of the I = m = 2 
mode, r is the distance to the source, M is the black hole mass, 
c is the speed of light, and G is the gravitational constant. 
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This is usually expressed in terms of the oscillation frequency 

fo — 'Si (0J22) /27r and the quality factor Q = 77/0/^ (1^22), 



h (t) = A 



GM 



-t/o*/Q cos (2nf t) . 



(2) 



Under the assumption that the waveform is completely known, 
we can implement the method of matched filtering (9), in 
which the data is correlated with a bank of signal templates 
parametrized by the ringdown frequency and quality factor. 
An analytic fit by Echeverria [ 1 1 to Leaver's numerical cal- 
culations ifTTl relates the waveform parameters to the black 
hole's physical parameters, mass M, and dimensionless spin 
factor, defined in terms of the spin angular momentum J, for 
the fundamental mode a = Jc/GM 2 : 




FIG. 1 : Venn diagram of the coincident detector times in hours. 
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where g(a) = 1 — 0.63(1 — a) 3 ^ 10 . Thus, if we detect the 
I = m = 2 mode, these formulae will provide the mass and 
spin of the black hole 0281 . 
The amplitude is given by 1 5 1 



A = 



eQ-*F(Q)-*g(a)- 



(5) 



The LIGO detectors are sensitive to gravitational waves in 
the frequency band of ^50 Hz to 2 kHz. This corresponds to 
a mass range of 1 1 to 440 Mq for a black hole with a = 0.9 
oscillating in its fundamental mode. Using a typical instan- 
tiation of the S4 noise power spectrum we can estimate the 
horizon distance, D^, the distance out to which a specified 
source with optimal location and orientation produces an SNR 
of 8 in the detector. We consider a spinning black hole with 
a = 0.9 and e = 1%. This quantity is shown as a function 
of mass, for each of the three LIGO detectors and the LIGO 
design sensitivity in Fig. [2] 



where F(Q) = 1 + 24 g 2 . In addition to a frequency and qual- 
ity factor dependence, the amplitude of the waveform also de- 
pends on the fraction of the final black hole's mass radiated as 
gravitational waves, e. This quantity scales with the square of 
the symmetric mass ratio 77, where r] = m\rri2 / (rrii + m2) 2 , 
and thus is largest for an equal mass binary lfT2l [131 . Nu- 
merical simulations of the merger of equal mass binaries have 
shown that approximately 1% of the final black hole's mass is 
emitted in gravitational waves |8|. In this search we do not at- 
tempt to evaluate e; we use the output of the filter to calculate 
the effective distance to a source emitting 1% of its mass as 
gravitational waves. The effective distance is the distance to 
an optimally located and oriented source. 



III. DATA SET 
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This search uses data from the 4 LIGO science run (S4), 
which took place between February 22 nd and March 24 th , 
2005. This yielded a total of 567.4 hours of analyzable 
data from the 4km interferometer in Hanford, WA (HI), 
571.3 hours from the 2 km interferometer in Hanford, WA 
(H2), and 514.7 hours from the 4 km interferometer in Liv- 
ingston, LA (LI). In this analysis we require that data be 
available from at least two detectors at any given time. This 
results in approximately 364 hours of triple coincidence and 
210 hours of double coincidence, as shown in Fig.[T] During 
S4, the LIGO detectors operated significantly below their de- 
sign sensitivity; this was attained in the subsequent science 
run, S5 fl4l. 



FIG. 2: The horizon distance versus mass and frequency for a black 
hole with spin of a = 0.9 and e = 0.01. From top to bottom, 
the curves show the horizon distance for the Initial LIGO reference 
design (black), the Hanford 4 km detector HI (red), the Livingston 4 
km detector LI (green), and the Hanford 2 km detector H2 (blue). 



During each science run there are times when disturbances 
couple into the data stream, introducing noise transients that 
can trigger a matched filter with high SNR. A careful study 
of correlations between LIGO data and auxiliary channels al- 
lowed us to identify data segments with an excessive noise 
transient rate, or with known artifacts. We refer to these as 
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data quality flags. These are then categorized according to the 
severity of the disturbance [ 15 1; triggers occurring during cat- 
egory 1 times are not analyzed as the excess of noise is likely 
to contaminate the estimation of the power spectral density. 
Data occurring during category 2 or 3 times are less prob- 
lematic, and thus to avoid excessive segmentation of the data, 
these are vetoed during the analysis. Any gravitational wave 
candidates occurring during category 4 times are cautiously 
examined. 



IV. PIPELINE 

When the signal is known, the optimal method of extract- 
ing the signal from Gaussian noise is matched filtering. LIGO 
data is non-Gaussian and while the method is still appropriate 
it is not sufficient to discriminate between signal and back- 
ground. We employ an analysis pipeline which was designed 
for this purpose, and closely resembles that of the inspiral 
searches described in lfl6l[T7l . Here we summarize the main 
steps. The first stage of the pipeline involves reading in and 
conditioning the data s(t) from each of the three LIGO de- 
tectors. We read in uncalibrated data, the differential arm 
length signal, at a sample rate of 16384 s -1 and down-sample 
to 8192 s _1 . This is converted to strain by applying the detec- 
tor response function. The one-sided power spectral density 
Sh(f) is calculated for each 2176 s long segment of the cali- 
brated data. The data is then broken further into sets of sixteen 
overlapping blocks, 256 s in length, and filtered using a bank 
of ringdown templates. The templates are positioned in f and 
Q according to the metric 1 1 8 1 



where 



ds 2 



1 



3+ 16Q 4 



Q 2 (l 



4Q 2 ) 2 
3 + 8( 



dQ 2 



3 + 4Q 2 



JoQ(l+4Q 2 ) 



d/ 2 



dQd/o 



(6) 



such that the maximum mismatch, ds 2 , between any point 
within the template bank and the nearest template is 3%. We 
search within the most sensitive portion of the LIGO fre- 
quency band, 50 Hz to 2 kHz and in quality factor between 
2 and 20. With these parameters we obtain a bank of 583 
templates with five different values of Q. The same tem- 
plate bank is used throughout the run. Numerical simula- 
tions |[T9l l20l |2T1 12211 have demonstrated that the maximum 
spin attained by the final black hole in a binary black hole 
merger is less than 0.96, corresponding to a quality factor 
of 8.5. In this search, we chose to cover a larger parameter 
space, and explore spin values between (Q — 2) and 0.994 
(Q = 20). Each filter has the form 



h(t) = cos (27r/ t) e^ ht ^ Q , 



< t < U 



(7) 



with a length of ten e-folding times, i max = 10 t, where r — 
Q/irfo- Filtering the data gives a signal to noise ratio (SNR) 



p{h) 



TP) 1 



(8) 



,h) 



S(f)h*(f) 
S h (\f\) 



df, 



(9) 



Here, the noise spectral density Sh(f) is the one appropriate 
for the data segment in question. For each filter, only trig- 
gers exceeding a pre-defined threshold, = 5.5 are retained. 
These are then clustered using a peak finding algorithm over a 
minimum of 1 s, with the loudest triggers for each filter writ- 
ten out to a file. Approximately 10 6 triggers were written out 
for each detector for the entire S4 run. 

In order to claim a detection of a gravitational wave ring- 
down we require coincidence between at least two detectors 
in the time of arrival of the signal and in the waveform param- 
eters. The time requirement is enforced first, only allowing 
triggers that appear within 4 ms of each other for co-located 
Hanford detectors and 14 ms for the Livingston-Hanford pairs 
(the two observatories are separated by 10 ms of light travel 
time) through to the next stage. The second coincidence test 
is based on the metric used to lay out the template bank, equa- 
tion |6]), with the size of the coincidence window around each 
template varying depending on its position within the parame- 
ter space. At this stage we also veto triggers occurring during 
times when category 2 and 3 data quality flags were on, and 
implement an amplitude consistency test between triggers ap- 
pearing in both Hanford detectors, such that the pair of trig- 
gers is retained only if the ratio of the HI effective distance to 
the H2 effective distance is less than 2. This results in list of 
coincident triggers found in two or three detectors, hereafter 
referred to as doubles and triples respectively. The final step is 
to cluster this coincidences over a time window of 10 s, retain- 
ing only the coincidences with the highest value of a detection 
statistic pus, a measure of the relative significance of the co- 
incidences. For triple coincidences (H1H2L1) the detection 
statistic was 



Pds = Pm + Pm + Pli ■ 



(10) 



The high level of false alarms associated with two detector 
coincidences (H1L1, H1H2, or H2L1), shown in Fig. [3] war- 
ranted a different detection statistic, 

p DS = min{p ifol + p iio2 , a x p ifol + b, a x p ifo2 + b}, (11) 

where suitable values of a and b were found to be 2 and 2.2 
respectively from an evaluation of the false alarm rate. (The 
evaluation of false alarm rates is discussed in section [V]) Note 
that in this ranking we do not take the square root of equation 
( |T0] > so as to emphasize the triple coincidences over the dou- 
bles. The ten coincidences with the highest value of the detec- 
tion statistic in each clustered coincidence category (triples or 
doubles) are followed up upon as detection candidates. 



V. TUNING THE SEARCH 

As with previous matched filtering searches in LIGO we 
implement a "blind analysis" of the data, that is, the con- 
straints are decided upon prior to looking at the full data set. 
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FIG. 3: Plot of background events (black dot) and simulated signal 
(red x's) found in double coincidence. The green lines denote the 
lines of constant detection statistic for two detector coincidences. 



FIG. 4: Hanford effective distance versus ringdown frequency /o for 
simulated signals found in coincidence. Simulated signals found in 
triple coincidence are marked as green x's, simulated signals found in 
double coincidence are shown as red dots and those signals found in 
double coincidence because the third detector was vetoed are circled 
in black. 



We chose values of the constraints that maximize the num- 
ber of ringdown signals recovered while minimizing the false 
alarm rate. The constraints under consideration include the 
SNR threshold, the size of the coincidence windows, the de- 
tection statistic, and the amplitude consistency test. 

The signal is modeled by adding simulated signals in soft- 
ware to the data stream and running the pipeline in the same 
manner as one would with the real data. Fig. [4] displays a plot 
of Hanford effective distance (the distance to an optimally 
located and oriented source) as a function of frequency for 
simulated signals that were found in all three detectors and 
in all combinations of two detectors (H1H2, HI LI, H2L1). 
The plot shows that the search is most sensitive to ringdown 
signals occurring in the 70Hz-140Hz band, where detector 
noise is lowest. 

The background, or false alarm rate, is estimated by shifting 
two data streams in time with respect to one-another and run- 



ning the pipeline described in section IV The LI data stream 



is slid fifty times in multiples of ±10 s and H2 is slid fifty 
times in multiples of ±5 s. As the time steps are much larger 
than the longest possible delay between detectors in receiv- 
ing a real gravitational wave signal, any coincidences found 
cannot be due to gravitational waves, and are therefore con- 
sidered false coincidences. We use these as estimates of the 
background due to false in- time coincidence. 

As a final sanity check before the search is unblinded we 
look at approximately one tenth of the data to ensure that the 
false alarm rate is consistent with our background studies. In 
order to avoid any potential bias in the tuning procedure af- 
fecting the upper limits, these data are excluded from the up- 
per limit calculation. 



VI. RESULTS 

Once the tuning is complete and all thresholds and cuts are 
finalized the pipeline described above is run on the full data 
set. Unblinding the search with the cuts and thresholds de- 
scribed above revealed no triple coincident events. A num- 
ber of double coincidences were found to satisfy the con- 
straints; however, these events are consistent with background 
as shown in Fig. [5] The loudest candidate events were sub- 
jected to further investigation, and in each case there was in- 
sufficient evidence that the coincidence could not be attributed 
to noise in the detector. 

We calculate an upper limit on the rate of ringdowns for 
a given population of black holes using simulated signals to 
evaluate the efficiency of the search, e(r), defined as the frac- 
tion of simulated signals detected in the analysis, as a func- 
tion of physical distance. The simulated signals were uni- 
form in orientation and direction. Given the expected high 
false alarm rate in double coincidence from background stud- 
ies, we set an upper limit from times when all three detectors 
were recording data, a total of T = 0.0375 years. Fig.[4]shows 
that the efficiency is strongly dependent on the ringdown fre- 
quency /o, and thus for the purpose of setting an upper limit 
we restrict the calculation to the most sensitive frequency 
band, 70Hz-140Hz, corresponding to black hole masses in 
the range 85 M Q -390 M Q . We added 5701 simulated signals 
in the frequency range of 70 Hz-140 Hz and between 0.5 Mpc 
and 10 3 Mpc in distance to the data, over ten runs. Of these, 
3010 were recovered in triple coincidence. Fig.[6]displays the 
efficiency as a function of distance for the frequency band of 
interest. (Note that the efficiency is not equal to 1 at nearby 
distances because we apply vetoes and treat these as a loss of 
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FIG. 5: Cumulative histogram of coincident triggers (red circles) and 
an estimate of the background (black x's): H1L1 doubles in triple 
time. 



efficiency rather than a loss in live-time.) From the efficiency 
we can calculate the effective volume we are sensitive to, V^f 



Vrff = 47T / e(r)r dr 



(12) 



For this band V e s = 2.6 x 10 6 Mpc 3 . This corresponds to a 
typical distance to a source of s=s 85 Mpc. The 90% confidence 
upper limit on the rate is given by 



^90% 



2.303 
TVT 



(13) 



eff 



which, for the 70Hz-140Hz band is 2.4 x 10~ 5 yr^Mpc" 3 . 
Even though our knowledge of the population of intermedi- 
ate mass black holes is poor, as discussed in section [I] we 
do know that the formation of stars in general scales with the 
blue-light luminosity emitted by galaxies, and as it is expected 
that the rate of binary coalescence (including ringdown) fol- 
lows the rate of star formation, we can work under the as- 
sumption that the rate of ringdowns also scales with blue light 
luminosity [24|. We introduce the cumulative blue luminosity, 
Cx, observable within the range of the search, 



C L = PL V e 



eff- 



(14) 



Cl has units of L\q, where Lio = 10 10 Lb,® and Lb,q 
is the B-band solar luminosity (1 MWEG is equivalent to 
1.7 L w ), and p L = (1.98 ±0.16) x 10~ 2 L W Mpc" 3 is the 
average blue luminosity density |251 . The 90% confidence 
upper limit on the rate in these units is given by 



^90% 



2.303 



(15) 
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FIG. 6: The efficiency of detecting simulated signal injections in 
triple coincidence in the 70 Hz- 140 Hz band, fit to a sigmoid. 



Due to our lack of knowledge about the population of 
black holes we assign no error to the astrophysical source 
population. Similarly, we assign no error to the waveform: 
comparison with numerical relativity results has shown that 
exponentially-damped-sinusoid templates perform well at de- 
tecting the signal and characterizing the black hole param- 
eters 0. We limit ourselves to evaluating systematic er- 
rors associated with the experimental apparatus and analy- 
sis method. The only error that we associate with the for- 
mer is calibration of the data, and with the latter is with the 
limited number of Monte Carlo simulations used to evaluate 
the efficiency. Errors in the calibration can cause the SNR 
of a signal to be incorrectly quantified, thereby introducing 
inaccuracies in the distance. As the efficiency is a function 
of distance and frequency, care has to be taken to adjust the 
efficiency curve appropriately. The fractional uncertainty in 
amplitude was found from calibration studies ||26l to be 5%. 
This results in an error of 3.5 x 10 5 Mpc 3 in V e s. The sec- 
ond source of error is due to the limited number of simu- 
lated signals in our Monte Carlo (MC) simulations to eval- 
uate the efficiency. Assuming binomial errors we calculate 
the variance of the efficiency er^C' which corresponds to an 
error in V e ts of 2.3 x 10 5 Mpc 3 . These errors are summed in 
quadrature, and multiplied by 1.64 to give a 90% confidence 
interval of 6.9 x 10 5 Mpc 3 . To obtain an upper limit we ap- 
ply a downward excursion to the effective volume and obtain 



R 



90% 



3.2 x lO^yr- 1 Mpc~ 3 = 1.6 x lO^yr^L 



lr-i 

10 



which evaluates to 1.2 x 10 3 yr 1 L 



10 



VII. CONCLUSION 

We have performed the first search for gravitational waves 
signals from black hole ringdowns in LIGO data, and demon- 
strated with simulated signals that this pipeline is an effective 
method of detecting such signals in coincidence between all 
three LIGO detectors. The search encompassed black holes in 
the mass range of 10A/ Q -500M Q , the regime of intermedi- 
ate mass black holes, with spins ranging from to 0.994. No 
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gravitational-wave events were found, and an upper limit of 
i?g % = 1.6 x 1CP 3 yr -1 !^ was placed on the rate of ring- 
downs from black holes formed from binary mergers, in the 
mass range of 85 M — 390 M . 

A search for black hole ringdowns in data from the fifth 
LIGO science run is currently underway. This two year long 
science run was the first at LIGO design sensitivity. With the 
increase in sensitivity of the LIGO detectors between the two 
runs, the ringdown horizon distance of S5 is greater than that 
of S4, allowing access to a greater number of potential sources 
from which a detection could be made. A further increase 
in sensitivity will come with Advanced LIGO, allowing us to 
detect compact binary coalescence to cosmological distances, 
and the improved sensitivity at lower frequency will make us 
sensitive to black holes with masses up to 2000 M or higher. 
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